Towards the chemical tuning of entanglement in molecular nanomagnets 
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Antiferromagnetic spin rings represent prototypical realizations of highly correlated, low- 
dimensional systems. Here we theoretically show how the introduction of magnetic defects by 
controlled chemical substitutions results in a strong spatial modulation of spin-pair entanglement 
within each ring. Entanglement between local degrees of freedom (individual spins) and collective 
ones (total ring spins) are shown to coexist in exchange-coupled ring dimers, as can be deduced from 
general symmetry arguments. We verify the persistence of these features at finite temperatures, and 
discuss them in terms of experimentally accessible observables. 
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From a magnetic perspective, most molecular nano- 
magnets (MNs) can be essentially regarded as spin clus- 
ters with dominant exchange interaction As such, 
they represent prototypical examples of correlated, low- 
dimensional quantum systems. It is thus tempting to 
consider how the wide tunability of their physical prop- 
erties, enabled by chemical synthesis, can allow to con- 
trol quantum entanglement In particular, differ- 
ent forms of entanglement can be possibly fine-tuned by 
chemical processing of well defined molecular building 
blocks, such as the controlled substitution of single mag- 
netic ions or the growth of supramolecular bridges, both 
of which have been recently demonstrated in Cr-based 
wheels 0, [||. The former process results in the intro- 
duction of magnetic defects in otherwise homogeneous 
molecules @, thus affecting correlations between their 
constituent spins. The latter one can instead induce weak 
exchange couplings between two or more MNs Q, so as 
to entangle their total spins @. 

The effect of magnetic defects on entanglement has 
been widely investigated in infinite systems of 1/2 spins 
0. With respect to these, MNs present some relevant dif- 
ferences, that make them peculiar model systems: they 
typically consist of s > 1/2 spins, are finite systems, 
and are characterized by point symmetries rather than 
translational invariance. The role played by these fea- 
tures [9]. as well as by the specific magnetic defects and 
supramolecular structures mentioned above, represents a 
specific reason of interest for the study of quantum en- 
tanglement in MNs. Here we refer to a prototypical class 
of molecules, namely the heterometallic octahedral &7M 
wheels, obtained by replacing one of the Cr ions in the 
parent Crs molecule with different transition metals M 
(If! [Tl| . We theoretically show how the magnetic defects 
spatially modulate spin-pair entanglement, and identify a 
clear dependence of such effect on the impurity spin sm- 
We finally investigate the interplay between intra- and 
inter-molecular entanglement in dimers of &7M rings, 
whose compatibility can be deduced from general sym- 
metry arguments. All these features are fully captured by 
local observables, such as exchange energy of individual 
spin pairs and partial spin sums, that are now available 
in direct geometry inelastic neutron scattering |12t Il3| . 

The dominant term in the spin Hamiltonian of the 
&7M rings is the antiferromagnetic exchange between 



FIG. 1: (color online) (a) Negativity N and entanglement of 
formation Ef for nearest neighboring spins of the Crs ring, 
as a function of temperature. The dotted lines are the gaus- 
sian fits of the reported points, (b) Temperature dependence 
of the entanglement witnesses Ws (red) and Wh (blue), nor- 
malized to the respective threshold values: W$ h = 12 and 
Wj^ = -18 J. The shaded areas correspond to the tempera- 
ture ranges where the two EWs detect entanglement. 



nearest neighbors [14] : 
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H = s k ■ s k+1 + J' s N ■ (s 



JV-l 



+ (1) 



fc=l 



where s k <N — scr = 3/2, sn — sm, and N — 8. Addi- 
tional terms, accounting for local crystal field and dipo- 
lar interactions, are typically two orders of magnitude 
smaller, and can be safely neglected in the present con- 
text. The spin sm introduced by the chemical substi- 
tution determines the total spin S of the ring ground 
state. In particular, the spin values for the elements 
M=Zn, Cu, Ni, Cr, Fe, Mn, are s M = 0, 1/2, 1, 3/2, 2, 5/2, 
corresponding to S — 3/2,1,1/2,0,1/2,1, respectively 

Mm 

In order to investigate the entanglement properties of 
the Cr7M molecules, we compute entanglement witnesses 
and measures [3, Il5j . An entanglement witness (EW) 
is an observable whose expectation value can exceed a 
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FIG. 2: (color online) Negativity of the neighboring spin pairs 
in the Cr7M rings. The values refer to the ground state, with 
M = S, of the spin Hamiltonian H . The inset shows a pic- 
torial representation of the entanglement localization induced 
by the defect M: the shaded areas highlight the most entan- 
gled spin pairs. 



given threshold only in the presence of (a particular form 
of) entanglement. For a ring consisting of N exchange- 
coupled spins s (with J > 0), the Hamiltonian itself can 
be regarded as an EW, being ( H) > —JNs 2 for any 
fully separable density matrix p [lfl, [TtJ ■ The violation 
of such inequality implies the presence of entanglement 
within the system, and in particular between pairs of 
neighboring spins. We generalize the above criterion to 
the case of a ring with a spin Sjv =/= Sk<N and two different 
exchange constants (Eq. [!}: 



W H = (H) > -(N - 2)Js 2 Cx - 2J's Cr s M = W J 



— IXlth 
H 7 



(2) 



being J, J' > 0. The full separability of the density ma- 
trix also implies a lower bound for the variances of the 
total-spin projections 0: E a = x , y J( S *) ~ ( S <*) 2 } > Na. 
In rotationally invariant spin Hamiltonians, the above 
quantity can be identified with the mag netic susceptibil- 
ity, up to a multiplicative constant [191 ] . We modify the 
above criterion to account for the presence of the mag- 
netic defect (Eq. [T|): 



E 

:=x,y,z 



[(^)-(Sa) 2 ] > (iV-l) SCr + SM 



w s h . 



(3) 



We note that, unlike Wh, Ws depends on correlations 
between all spin pairs, and allows to detect forms of en- 
tanglement that don't show up in the two-spin density 
matrices. 

While EWs represent a practical means for the detec- 
tion of entanglement, its quantification requires entangle- 
ment measures [3J. In particular, entanglement between 
the spin pairs (si,Sj) is hereafter quantified by the neg- 
ativity J\f of their reduced density matrices py. This 
is defined as: N(pij) — (J2k \^k\ ~ ^fe)/2, where are 
the eigenvalues of the partially transposed density matrix 

PiJ' 

Homometallic ring — The parent molecule of the het- 
erometallic CtjM rings is Cr§ [2(| ■ Its energy spectrum is 



characterized by an S = ground state, separated from 
the first excited multiplet by a gap A ~ 0.559 J [2l[. As 
for the case of an 8 qubit ring [22], |23J , spin- pair entan- 
glement in the absence of an applied magnetic field is 
limited to nearest neighbors: we find in fact that Af = 
at all temperatures for any other spin pair. The tempera- 
ture dependence of Af(pi,i+i) presents instead a finite and 
nearly constant value for T/J < A, followed by a smooth 
decay at higher temperatures [Fig. [T](a), black squares]. 
The threshold value of the temperature at which M van- 
ishes is given by: T^f ~ 1.58J. The same qualitative be- 
havior is found for the entanglement of formation Ep 0] , 
here reported as a benchmark (red squares). The witness 
Whi like Af, only reflects entanglement between neigh- 
boring spins, and in fact provides a similar value for the 
threshold temperature: W H < W# for T < T H ~ 1.5J 
[Fig. Q] (b)]. The witness Ws reveals instead further 
forms of entanglement, that persist up to higher tem- 
peratures, being W s < Wf for T < T s ~ 2.75J. We 
finally note that values of Wh/J lower than -20.65 and 
-21.8 imply the presence in the system state of three- 
and five-spin entanglement, respectively [24T ]. In the Cr§ 
ring, these thresholds are exceeded for T/J lower than 
1.05 and 0.74. 

Heterometallic rings — The chemical substitution of 
one Cr with an M ion introduces a magnetic defect in 
the molecule. This consists in the replacement of an 
scr with an sm ^3/2 spin, and - in principle - in the 
modification of the exchange coupling between s% and its 
nearest neighbors (Eq. [T|). In the following we focus on 
the case J' = J, which is compatible with the current 
estimates [ll[, an( i allows to isolate the dependence on 
sm of spin-pair entanglement. In particular, we compute 
the negativity Af(pi~i.i) to quantify the ground-state en- 
tanglement, and the local witnesses Wat '■ 



W Hi = J(si-i ■ > -Js{ 



-lSj 



(4) 



Violation of the above inequality implies entanglement 
specifically between s^_i and s,. The witnesses Wh { , that 
can be regarded as the local versions of Wh = Yli=i Wat , 
are used hereafter to estimate the threshold tempera- 
tures up to which entanglement persists in the various 
molecules. 

The negativity corresponding to the ground state of 
H is reported in Fig. [2] The translational invariance 
that characterizes the case of Crg (black squares) is re- 
placed by a strong spatial modulation. The negativity of 
the spin pairs with i = 1, 3, 6, 8 is an increasing function 
of Sm, while for i = 2,7 the opposite occurs. For each 
molecule, the oscillations of A/^Pi-i^) as a function of 
i are paralleled by those of WHi ( n °t shown). In fact, 
they can be intuitively explained in terms of competi- 
tion between the non-commuting exchange operators of 
H, since the state that maximizes entanglement between 
Si and i coincides here with the one that minimizes 
their exchange energy (25| . In particular, large values of 
the impurity spin tend to minimize Wh x (and thus to 
maximize entanglement of sm with s\), at the expense of 
the competing term Wh 2 (and of entanglement between 
si and S2). This in turn favors exchange interactions and 
quantum correlations between S2 and S3 , and so on. The 
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FIG. 3: (color online) Threshold temperatures of the local 
witnesses WH t as a function of the impurity spin sm: (a) 
M-Cr pair (i — 1); (b) three inequivalent Cr-Cr pairs i — 
2, 3, 4 (blue, orange, green). The grey curve, reported in both 
panels, gives the threshold temperatures of the global witness 
Wh- (c) Threshold temperatures for the witnesses Ws n as a 
function of the spin number n and of the impurity spin sm- 



inverse applies to rings with sm < so- 

These general trends are confirmed by the finite- 
temperature results. In Fig. [3] (a,b) we report the tem- 
peratures below which the witnesses Wh ( violate the in- 
equalities Eq. |4l thus detecting entanglement in the cor- 
responding spin pair. The spin pairs that exhibit the 
strongest dependence on sm are the ones that include or 
are next to the magnetic defect: in particular, Th 1 in- 
creases linearly with sm [panel (a), red squares], while 
Tjj 2 decreases monotonically [panel (b), blue squares]. 
The other two inequivalent spin pairs also display a 
monotonic - though weaker - dependence on sm (green 
and orange squares). These different behaviors cannot 
be appreciated in Th(sm) (grey curve in both the upper 
panels). In fact, the threshold temperature of the global 
witness Wh increases monotonically for finite values of 
the impurity spin, and systematically underestimates the 
persistence of entanglement, being > Th for at least 
some i in all molecules. 

In order to bridge local (two-spin) and global witnesses, 
we introduce the EWs Ws n , corresponding to the vari- 
ances of partial spin sums: 



W Sn 



a=x,y,z 



[((s n , a y)-(s n , a y]>y si = wf. (5) 



1=1 



where S n = Ya=i Sl anc ^ tnus Ws n = Ws- Their thresh- 
old temperatures T$ n are reported in Fig. [3Jc) for even 
values of n. The most prominent feature is the signifi- 
cant increase of T$ n with n, for all the Cr-iM molecules. 
The same feature shows up within n— spin subensembles 
with odd n (not shown). Entanglement for Tjj < T < T5 
is thus progressively averaged away by the partial traces 
that are performed to derive the reduced density ma- 
trices p n = Tr, i+ i a{p} of decreasing spin number n. 
Based on calculations performed for the 8-qubit ring, we 
conjecture that the density matrix p(T) in such temper- 
ature range is given by the mixture of terms that include 
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FIG. 4: (color online) Negativity of neighboring spin pairs 
belonging to CryM ring dimers, with M=Ni (a) and M=Cu 
(c): AA/i = 1— A/i(a)/A/i(0). (b,d) Trace distance A between 
the reduced density matrices pf-\ : i{ot) and pt-i,i{0). The 
correspondence between spin pairs and colors is reported in 
panels (a,c). 



two-spin entanglement between different pairs. In other 
words, p(T) would be a mixed 2-producible state [26j j . 
whereas multi-spin entanglement is limited to T < J (see 
discussion on the Crg molecule). As to the dependence 
on the impurity spin, this is most significant for n = 2 
and n = 4. In particular, Ts 2 decreases monotonically 
with sm, like Th 2 (that refers to the same spin pair), 
while Ts s = T$ displays the opposite behavior. 

Ring dimers — The CryM rings are composite quan- 
tum systems, whose subsystems are represented by the 
constituents magnetic ions. However, they can also be re- 
garded as building blocks of supramolecular assemblies, 
consisting of molecule dimers or oligomers Q. Here, the 
intermolecular exchange typically mixes states belonging 
to the ground S multiplet of each ring, and can thus 
entangle collective degrees of freedom, such as the total 
spin projections of the MNs H . 

In the following, we investigate the interplay between 
individual- and collective-spin entanglement in a proto- 
typical state of a CryNi-ring dimer: 

\9%(a)) = a\l/2, -1/2) - (1 - a 2 ) 1/2 | - 1/2, 1/2), (6) 

where A and B label the two rings (with Sa/b — 1/2), 
and the components \Ma, Mb) are labelled after the val- 
ues of the ring-spin projections (Sf and Sf). Entan- 
glement between these collective spins can be quanti- 
fied by the negativity: Nab = J^Q^ab)) = a^/l-a 2 . 
The state I^ab) thus varies from a maximally entangled 
(A/as = 1/2) to a factorizable state {Nab — 0) as a 
ranges from l/y2 to 1. Entanglement between two in- 
dividual spins, sf and s^, of ring A is quantified by the 
negativity of the two-spin reduced density matrix pf, : 



a p. 



(1 - o?)p 



-1/2 

ij 



(7) 



Here, pf^ A denotes the two-spin density matrix corre- 



sponding to the ground state of the single-ring Hamilto- 
nian (Eq. [1]) with total-spin projection Ma- 
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Entanglement between the total spin projections of A 
and B (a < 1) results in a mixing of p\( 2 and p^ 1 ^ 2 (Eq. 
[7]) . This generally tends to reduce entanglement between 
Si an d Sj - due to the convexity of entanglement mea- 
sures [151 ] - making quantum correlation between individ- 
ual and collective degrees of freedom mutually exclusive. 
As shown in Fig. |4l such reduction is however very lim- 
ited in the case of exchange-coupled rings. In fact, as a 
varies from l/y/2 to 1, the relative change of the nega- 
tivity M{pi-\,i) is below 2% for all the spin pairs (panel 
a), in spite of the fact that the reduced density matri- 
ces pi—i t i vary significantly with ajb). Such variation is 
quantified by the trace distance [271 ] between pij(a) and 
the reference state Pij(l)- 

Ai (a) = (l/2)Tr 0^- (a) - Pij (l)]t [ Pij (a) - Pij (1)] . (8) 

Entanglement between individual spins within each 
CryNi ring (sf and s^), is thus fully compatible with 
that between collective spins (Sf and Sf) : such as that 
induced by weak intermolecular exchange. 

The same applies to different ring dimers, such as that 
formed by two CryCu rings (Sa = Sb = !)• Here, we 
consider the prototypical state 

= [(1 - « 2 )/2] 1/2 (|l, -1) + | - 1, 1)) - a\0, 0). 

(9) 

This passes from maximally entangled (Nab = 1) to fac- 
torized (Nab = 0), as a varies from l/-\/3 to 1, being 
Nab — cty / 2(l — a 2 ) + (l— a 2 )/2. In this same range, the 
negativity displays a very limited decrease for all pairs of 
neighboring spins (panel c), also for significant values of 
the trace distance (d). 

Symmetry arguments suggest that such compatibility 
between intra- and inter-molecular entanglement is more 



general. In fact, the entanglement witnesses Whi are 
scalar operators. As results from the Wigner-Eckart the- 
orem [281 ] . their expectation value is thus identical for 
all the states that form an irreducible representation of 
the rotational group, such as the eigenstates belonging to 
the ground S multiplet of the spin-ring Hamiltonian H 
(Eq. [I]). As a consequence, if any of the inequalities Eqs. 
I2l4l is violated by a single-ring ground state, it's also vi- 
olated by ring-dimer singlet states such as \^^ B (1/V2)} 
and ^^(l/V^i)), where and Sf are maximally en- 
tangled and the single-ring density matrix is a mixture 
of all the ground multiplet states. The same argument 
applies to the witnesses Ws n , provided that their state is 
rotationally invariant: in this case, (S n ,a) = 0, and Ws n 
can be identified with the expectation value of the scalar 
operator S 2 . 

In conclusion, we have shown how the introduction of 
magnetic defects in the heterometallic Cr-^M rings intro- 
duces a strong spatial modulation of pairwise entangle- 
ment, that persists at finite temperatures. This suggests 
that suitably combined chemical substitutions can rep- 
resent an effective means for engineering entanglement 
in molecular systems. Besides, we quantitatively show 
that in ring dimers entanglement between individual and 
collective spins can coexist, and deduce the generality of 
such property from symmetry arguments. The discussed 
features are fully captured by local observables acting as 
entanglement witnesses, such as the exchange energy of 
spin pairs and the variances of partial spin sums, that are 
accessible to direct geometry inelastic neutron scattering. 
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